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Abstract 

Global warming can have substantial impacts on the phenological and growth patterns of alpine and Arctic species, 
resulting in shifts in plant community composition and ecosystem dynamics. We evaluated the effects of a six-year 
experimental soil warming treatment (+4°C, 2007-2012) on the phenology and growth of three co-dominant dwarf shrub 
species growing in the understory of Larix decidua and Pinus uncinata at treeline in the Swiss Alps. We monitored vegetative 
and reproductive phenology of Vaccinium myrtillus, Vaccinium gaultherioides and Empetrum hermaphroditum throughout 
the early growing season of 2012 and, following a major harvest at peak season, we measured the biomass of above-ground 
ramet fractions. For all six years of soil warming we measured annual shoot growth of the three species and analyzed ramet 
age and xylem ring width of V. myrtillus. Our results show that phenology of the three species was more influenced by 
snowmelt timing, and also by plot tree species {Larix or Pinus) in the case of V. myrtillus, than by soil warming. However, the 
warming treatment led to increased V. myrtillus total above-ground ramet biomass (+36% in 2012), especially new shoot 
biomass (+63% in 2012), as well as increased new shoot increment length and xylem ring width (+22% and +41%, 
respectively; average for 2007-2012). These results indicate enhanced overall growth of V. myrtillus under soil warming that 
was sustained over six years and was not caused by an extended growing period in early summer. In contrast, £ 
hermaphroditum only showed a positive shoot growth response to warming in 201 1 (+21%), and V. gaultherioides showed 
no significant growth response. Our results indicate that V. myrtillus might have a competitive advantage over the less 
responsive co-occurring dwarf shrub species under future global warming. 
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Introduction 

Alpine and Arctic ecosystems are predicted to be especially 
vulnerable to global warming [1,2] because plant growth and 
performance in these environments are strongly constrained by 
low temperature, short growing seasons and frequent freezing 
events during the snow-free period. Projected increases of 1 .8— 4°C 
in the global mean surface air temperature by the end of the 2 1 st 
century [3] can cause dramatic changes in community composi- 
tion since co-occurring species respond differently to climate 
variations [4-6]. Community changes are of particular importance 
in alpine and Arctic ecosystems due to their potential effects on 
climate feedbacks through shifts in plant community composition 
and plant species cover [7,8]. 

Experimental studies have shown that responses of plants 
growing at high latitude and high elevation to warming are 
species-specific [9-11], demonstrating that understanding how 
plant community dynamics might change with higher tempera- 



tures can only be achieved through assessing responses of 
individual species. For some species and study sites, positive 
growth responses to warming were transient, stopping some years 
after the experiment started [12,13]. In contrast, other studies 
found a lag in plant growth response to experimental warming or 
inconsistent responses over time [6,14]. Therefore, species-specific 
studies lasting several years are crucial for understanding 
community responses to global warming. 

Phenological timing is an important factor for plant fitness and 
abundance, and it can be highly responsive to temperature ([15] 
and references therein). For instance, advanced leaf phenology in 
response to warming could lead to an extension of the 
photosynthetically active season, which, in turn, could lead to 
greater carbon gain and ultimately growth for plants [16,17]. 
Advanced phenology of high-latitude and high-elevation species 
has occurred with experimentally increased growing season 
temperatures in numerous studies where snowmelt date was not 
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modified [10,18-22] (but see [23]). In a meta-analysis of alpine 
and Arctic studies, Arft et al. [12] suggested a relationship between 
advanced leaf bud burst under warming and increased vegetative 
growth in alpine and Arctic species. However, most studies have 
focused on phenology and growth separately, whereas the 
relationship between phenology and growth in response to 
warming at an individual/ ramet level has rarely been reported. 

Shrubs are an essential component of alpine and Arctic 
ecosystems, and several studies have reported an increase in their 
cover due to climate warming [6,14,24-27]. This "shrubification" 
strongly impacts ecosystem dynamics because shrubs modify 
patterns of snow accumulation, decrease albedo and modify 
carbon storage and nutrient cycling by greater biomass accumu- 
lation and more recalcitrant litter (see review by Myers-Smith et al. 
[8])- 

Shrub growth responses to global warming in alpine and Arctic 
ecosystems have been assessed in many studies through long-term 
observations [5,28], experimental manipulations [6] and/or 
studies using gradients across elevation or latitude [29,30]. Most 
previous warming experiments have measured above-ground 
primary production by means of new shoot growth or standing 
biomass [11,21,31-33]. In contrast, relatively few studies have 
focused on shrub secondary growth [13,34,35] since acquiring 
these data requires destructive techniques. Therefore, it remains 
unclear if the reported increases in biomass production or in shoot 
increment length are caused by shifts in biomass allocation or if 
greater overall biomass production occurs. 

The Stillberg treeline research area in the Swiss Central Alps 
hosted a six-year soil warming experiment (2007-2012). A short- 
term investigation of the dominant ericaceous dwarf shrub species 
growing in the understory of treeline trees after three years of 
warming showed that Vaccinium myrtillus responded to higher 
temperatures with increased shoot increment length but that 
Vaccinium gaultherioides and Empetrum hermaphroditism did not show a 
growth response [1 1]. Moreover, there were only few indications 
for changes in early-season vegetative phenology in these species 
under warming [36,37]. In the study presented here, we 
conducted a detailed investigation of dwarf shrub phenological 
and growth responses during the summer of 2012, the final (sixth) 
year of the soil warming experiment. A final destructive harvest 
conducted at the peak of the vegetation period allowed us to 
investigate warming effects on secondary growth and biomass 
allocation for the first time. Our aims were (i) to determine 
whether a positive growth response in Vaccinium myrtillus was 
maintained after six years of soil warming, and whether Vaccinium 
gaultherioides and Empetrum hermaphroditum showed any delayed 
increases in growth after this full experimental period. By tracking 
the vegetative and flowering phenology of these dwarf shrubs, we 
also aimed (ii) to relate possible lagged changes in their phenology 
to changes in growth. Finally, we analyzed the widths of V. myrtillus 
growth rings with the aim of (iii) assessing whether the increases in 
shoot increment length under soil warming previously reported for 
this species represented an increase in the entire above-ground 
biomass of the ramets or if there was simply a shift in biomass 
allocation. 

Material and Methods 

Study site 

The study site was located within the Stillberg treeline site in the 
Central Alps (Davos, Switzerland, 9° 52'E, 46° 46' N). Stillberg is 
a 5-ha long-term afforestation research area where tree seedlings 
were planted into the intact dwarf shrub community in 1975 by 
the Swiss Federal Institute for Forest, Snow and Landscape 



Research (WSL). Climate data measured by a WSL meteorolog- 
ical station located within the research area (2090 m a.s.l.) 
indicated a mean annual precipitation of 1 1 55 mm and mean 
annual air temperature of 2.1°C from 1975 to 2012. For the same 
period, the main growing season months (June-August) had a 
mean precipitation of 444 mm and a mean air temperature of 
9.2°C. See Table SI for details on climate data over the study 
years (2007-2012). 

No specific permits were required for this location and activities 
and the field studies did not involve endangered or protected 
species. 

Experimental design 

The experiment consisted of 40 hexagonal 1.1 m 2 plots, 20 with 
one Pinus mugo ssp. uncinata (DC.) Domin individual in the center 
and 20 with one Larix decidua Mill, individual in the center. The 
plots were situated within an area of 2500 m 2 on a NE-exposed 
25-30° steep slope at 2180 m a.s.l. at or slighdy above the current 
treeline in the region [38,39]. The trees were sparsely distributed 
without forming a closed canopy; therefore, each plot contained a 
single tree surrounded by a dense cover of understory vegetation 
consisting mainly of the co-dominant dwarf shrub species Vaccinium 
myrtillus L., Vaccinium gaultherioides Bigelow (group V. uliginosum agg.) 
and Empetrum nigrum subsp. hermaphroditum (Hagerup) Bocher 
(referred to hereafter as Empetrum hermaphroditum). Further details 
about understory species composition were reported by Martin et 
al. [36] and Dawes et al. [11]. 

A free air CO z enrichment (FACE) experiment was started after 
snowmelt in early June 2001 and applied during each snow-free 
season for nine years (ending in 2009). The 40 plots were assigned 
to 10 groups of four neighboring plots (two Pinus and two Larix 
trees per group) and half of these groups were randomly assigned 
to an elevated C0 2 treatment while the other half served as 
controls (see a detailed description of the setup and performance of 
the CO z enrichment facility in Hattenschwiler et al. [38] and in 
Dawes et al. [40]). Dwarf shrub responses to the C0 2 enrichment 
have been reported in detail by Dawes et al. [11,41]. A soil 
warming treatment was added to the experiment in spring 2007 
and was applied during each snow-free season until early August 
2012. Within each of the 10 C0 2 treatment groups, one plot of 
each tree species was randomly selected and assigned a soil 
warming treatment, either control or warmed, and the second plot 
was assigned the other treatment, yielding a balanced split-split- 
plot design with a replication of five individual plots for each 
combination of C0 2 level, soil warming treatment and tree 
species. Therefore, from 2007 to 2009 the experiment included 
both C0 2 enrichment and soil warming, whereas from 2010 to 
2012 the treatment consisted of soil warming only. The soil 
warming treatment was applied using 420 W-heating cables laid 
on the ground surface underneath the dwarf shrub layer, with a 
5 cm distance between neighboring cables [42]. Each year, the 
heating was turned on immediately after snowmelt and turned off 
just before the site was covered in snow for the winter to avoid an 
interaction between soil temperature and snow cover duration. 
The warming treatment increased the growing season mean soil 
temperatures at 5 cm depth by 3. 1 to 4.4°C over the 6 seasons of 
heating, with a difference of +3.5°C in 2012 (HOBO U23 Pro v2 
dataloggers, Onset Computer Corporation, Bourne, MA, USA). 
Air temperature was warmed within the dwarf shrub canopy 
(0.9°C at 20 cm above ground [42]). The warming treatment had 
a slight drying effect on the soil organic layer (details in Hagedorn 
et al. [42] and Dawes et al. [11]). However, the soil matric water 
potential at 5 cm depth was always above -300 hPa in all plots, 
indicating overall very moist soil conditions [11,43]. 
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Snowmelt date of each plot from 2007 to 201 1 was defined as 
the date in spring when soil temperatures at 5 cm depth rose 
sharply from values near 0°C, supported by visual estimations in 
the field and photographs. In 2012, snowmelt date was determined 
by visual estimations as the date when 50% of the plot was snow 
free and the ramets selected for detailed investigation were 
uncovered. This way, the estimation was more accurate for the 
study of the specific ramets selected. 

Phenology 

At the time of snowmelt in 2012, we marked five ramets of V. 
myrtillus, V. gaultherioides and E. hermaphroditum in each plot, 
excluding the area within 10 cm of plot borders to avoid potential 
edge effects. All plant measurements were made on these ramets. 
Vaccinium myrtillus was present and abundant in every plot, V. 
gaultherioides was present in 35 of the 40 plots, and E. hermaphroditum 
was present in 26 of the 40 plots, which was sufficient replication 
to assess treatment effects. 

We monitored the phenology of the three study species between 
the start of snowmelt (day of year 142) and peak growing season 
(day of year 212). We visited all plots and monitored all ramets 
every 2 to 4 days at early growing season and every 4 to 6 days 
after flowering. For each marked ramet, we recorded the date 
when it entered the following phenophases (some of which 
occurred at the same time): (1) burst of first vegetative bud, (2) 
first leaf starting to unfold, (3) first leaf fully expanded, (4) start of 
shoot elongation, (5) burst of first flower bud, (6) anthesis, and (7) 
first flower withered. In early August, we harvested all vegetation 
from the experimental plots. Previous phenological data in the 
same plots from an entire snow-free season [44] confirm our 
observation that dwarf shrub vegetative development was com- 
pleted at the time of the harvest, whereas fruits were not mature 
and leaf senescence had not started. 

Shoot increment length and above-ground biomass 

To track shoot growth responses during the whole soil warming 
experimental period, shoot increment length of the three species 
was measured from 2005 to 201 1 on the longest branch of five to 
seven ramets in the field (2005-2009 data presented in Dawes et 
al. [11]). During early August 2012, we harvested all the marked 
ramets in every plot, clipping them at ground level. For E. 
hermaphroditum, we measured the new shoot increment length on 
the longest branch for every ramet before detaching and drying 
them for biomass measures. For all three species, we separated 
new shoots and leaves from the rest of the ramet while they were 
still fresh and dried them at 60°C for 48 h. For the two deciduous 
Vaccinium species, we weighed new shoots and leaves separately 
and measured the length of three new shoots per ramet, whereas 
for E. hermaphroditum, we weighed the small needle-like leaves and 
new shoots together. For V. myrtillus and V. gaultherioides, we 
counted the number of new shoots per ramet and calculated the 
average mass per individual shoot. For E. hermaphroditum, we 
counted the number of new shoots from a subsample of ramets 
and used the mass-count relationship to estimate the number of 
branches for the remaining ramets in each plot. For all three 
species, we dried the remaining ramet material (designated as 
"main stems", > 1 year old) at 60°C for 24 h and weighed it. 

Growth rings of Vaccinium myrtillus 

We made cross-sections of 20 um thickness from the basal 
1.5 cm of dried V. myrtillus stems using a sledge microtome (WSL- 
Lab-microtome). We stained sections with a mixture of Safranin 
and Astrablue to emphasize the growth ring structure. For 
dehydrating the sections for preservation, we rinsed them with 



increasingly concentrated ethanol solutions (75%, 96%, 100%), 
immersed them in Xylol, imbedded them in Canada-Balsam and 
dried them at 60°C for 24 hours [45] before photographing them 
at x20-x200 magnification through a microscope with a digital 
camera (Canon EOS 650D on Olympus BX41 microscope; Fig. 
SI). We used the images to visually count rings and to measure 
xylem ring widths in four radii per section with the program 
WinCELL [46]. We excluded 23 out of 200 ramets from the 
statistical analyses because wood was damaged or broken or rings 
were not distinguishable. Ramets were not old enough for rings to 
be analyzed statistically using specialized dendrochronological 
software, so we visually cross-dated the samples to find possible 
missing rings. Latewood had not formed in the 2012 growth rings, 
indicating that secondary growth was not finished when we 
harvested the ramets. 

Statistical analysis 

We applied linear mixed effects models fitted with the restricted 
maximum likelihood estimation method (REML) to assess 
treatment effects on phenology and growth parameters. We used 
likelihood ratio tests to determine whether the previous C0 2 
treatment and interactions between C0 2 and the other treatments 
contributed significandy to the model fits as fixed effects. Previous 
investigations already reported a lack of soil warming x CO2 
interactive effects [1 1,41] and we did not find significant persistent 
effects of the C0 2 enrichment; therefore, we pooled across C0 2 
treatments for all final analyses. Statistical models for phenology 
and biomass of the study species included warming treatment, plot 
tree species and their interaction as fixed effects. We included 
snowmelt date as a covariate when it contributed significandy to 
the model fits, i.e. in tests of the effect of soil warming and plot tree 
species on phenology. 

Shoot increment length and growth ring width were analyzed as 
repeated measures and included treatment year (categorical 
variable) and all interactions between year, warming and plot 
tree species as additional fixed effects (after Dawes et al. [1 1]). 
Measurements averaged over 2005 and 2006 were included as a 
covariate to account for differences in shoot increment length or 
ring width before the warming treatment started. When a 
significant interaction between soil warming and year was found, 
we additionally tested the effect of soil warming with separate 
analyses for each individual year. We accounted for violation of 
independence of residuals from different treatment years by 
implementing the residual autocorrelation structure corARl [47]. 
As growth rings were not completely finished in 2012, we excluded 
this year from the repeated measures analysis of growth ring 
widths. We tested for Pearson's correlations between shoot 
increment length and above-ground biomass, as well as between 
shoot increment length and xylem ring width in 2012. We did not 
extend this analysis to previous years because ramets harvested in 
2012 were not the same as those randomly selected for shoot 
increment measurements in previous years. 

For all statistical analyses, the random effects structure reflected 
the hierarchy of the split-split-plot experimental design, with 
measurements made on ramets in 40 individual plots, nested 
within 20 soil warming treatment groups, nested within 10 C0 2 
treatment groups. For all analyses, we visually checked assump- 
tions of normality and homoscedasticity of the residuals. We log- 
transformed response variables when necessary to reach these 
assumptions. Effects were considered significant at P<0.05 and, 
because of the relatively low replication, we considered P&0.05 
but <0. 10 as marginally significant. We performed all the analyses 
with R version 2.15.2 [48] using the nlme package [47]. 
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Results 

Abiotic conditions 

In 2012, snow melted from the experimental plots between 21 
and 31 May (days of year 142 and 152). There were no differences 
in the date of snowmelt between soil warming treatments 
(F J>9 = 0.64, P= 0.447) or plot tree species (F lyls = 1.04, 
P= 0.32 1). From 2 1 May to 3 August, when plants were harvested, 
the mean air temperature at the meteorological station was 9.7°C 
and temperatures ranged from -0.4 on 12 June to 21.9°C on 30 
June. 

Phenology 

The soil warming treatment did not significandy affect the 
vegetative phenology of the three dwarf shrub species (Fig. 1). 
However, we found marginally significant warming effects on their 
flowering phenology. Vaccinium myrtillus and V. gaultherioides showed 
a slightly earlier flower anthesis (by 1 and 2 days, respectively) in 
warmed plots than in unwarmed plots (V. myrtillus: F ]9 = 4:.3, 
P = 0.067, Fig.lA; V. gaultherioides: i^ j7 =4.41, P= 0.074, Fig. IB). 
Vaccinium myrtillus flowers also withered slighdy earlier with 
increased soil temperatures (by 2 days, -F/ j g=4.03, P= 0.076). 

The tree species present in the plot had a greater effect than soil 
warming on V. myrtillus vegetative and reproductive phenology 
(Fig. 1A). The start of leaf unfolding, start of shoot elongation, 
flower bud break and flower anthesis took place 1 day earlier 
under larch than under pine (P<0.011), and leaf full expansion 
occurred 3 days earlier (Fj jy =29.88, P<0.001). There was a 
marginally significant effect of the tree species x warming 
interaction on the start of shoot elongation for V. myrtillus, which 
was slighdy earlier (1 day) in warmed plots under larch than in 
unwarmed plots under pine (F 117 = ?>.7 , P= 0.071). The phenol- 
ogy of V. gaultherioides and E. hermaphroditum did not show significant 
differences between plots with the two different tree species 
(Fig. 1B,C). 

Most of the phenological phases for V. myrtillus and V. 
gaultherioides occurred earlier in plots with an earlier snowmelt 
date (P<0.07). However, V. myrtillus anthesis and V. gaultherioides 
flower bud break and withering did not show a relationship with 
snowmelt timing. Empetrum hermaphroditum flowering phenophases 
and start of shoot elongation were also related to snowmelt date 
(P<0.04), but vegetative bud break was not. 

Shoot increment length 

Soil warming had a significant positive effect on the annual 
shoot increment length of V. myrtillus (mean increase of 22% over 
2007-2012, F 1:9 = 41.38, P<0.001, Fig. 2A). Treatment year also 
had a significant effect (F 5>180 = 3.67, P— 0.004) but the warming x 
year effect was not significant (i^iso = 1-89, P= 0. 100). Neverthe- 
less, the response was larger in the two last study years (shoot 
increment length 34% greater in warmed plots than unwarmed 
plots for 2011-2012) than in the previous years (18% greater in 
warmed plots than unwarmed plots, averaged over 2008-2010). 

The annual shoot increment length of V. gaultherioides was not 
significandy affected by soil warming (F 1>9 = 0.99, P= 0.346, 
Fig. 2B) although there was a trend of increase with warming in 
the last year. Treatment year had a significant effect on V. 
gaultherioides shoot increment length {F s ^i = 6.93, P<0.001) but 
the warming x year interaction was not significant (F s j 7] = 0.5, 
P= 0.776). Empetrum hermaphroditum annual shoot increment length 
showed a significant effect of treatment year [F 5 j 2 n- 25.97, P< 
0.001, Fig. 2C) and the plot tree species x year interaction 
(F5J2X- 2.55, P= 0.031), and also a marginally significant effect of 
warming (F; 9 =4.83, P= 0.056) and the warming x year 



interaction (F 5>12 s = 2.21, P= 0.057). Analyses of individual years 
showed that E. hermaphroditum only had a significant (positive) shoot 
growth response to soil warming in 20 1 1 (+2 1 % increase in 
warmed plots compared to unwarmed plots, F 1 ^= 10.01, 
P— 0.016) and that the shoot increment length of this species 
was greater under pine than under larch in 2008 (marginally 
significant, F U1 = \S)\, P= 0.070) and 2009 {F ul) =5.Q7, 
P= 0.048). 

Shoot increment length averaged over 2005 and 2006 (pre- 
treatment covariate) positively influenced the length of new V. 
gaultheroides shoots during the 2007-2012 warming period 
{Fjj7= 18.61, P<0.001), and to a lesser extent that of V. myrtillus 
and E. hermaphroditum (both marginally significant, i ? ; > ;7=3.28, 
P= 0.088 and i^„ = 4.58, P= 0.052, respectively). In 2012, new 
shoot increment length showed a positive correlation with new 
shoot biomass for each of the three species (V. myrtillus Pc =0.53, 
P<0.001; V. gaultherioides R 2 = 0.'i9, P<0.001; E. hermaphroditum 
R 2 = Q.27, P<0.001). 

Above-ground biomass 

At the time of harvest in August 2012, V. myrtillus ramets showed 
a 54% greater total leaf biomass {F 1>9 = 15.52, P= 0.003), a 63% 
greater total new shoot biomass (F 19 = 14.15, P= 0.005) and a 
26% greater main stem (stems >1 year old) biomass (F lt g — 7 .82, 
P— 0.021) in warmed plots than in unwarmed plots (Fig. 3A). 
Moreover, the main stem biomass of V. myrtillus was 35 % greater 
under pine than under larch (Fj j 8 = 5.15, P= 0.036). Although 
differences in ramet above-ground biomass between warming 
treatments were, on average, larger in plots with larch than in plots 
with pine, there were no significant tree species x warming 
interactions for any of the plant parts (P>0.77). Vaccinium 
gaultherioides and E. hermaphroditum did not show significant 
differences among treatments in the above-ground biomass for 
any of the plant parts analyzed (Fig. 3B,C). The three dwarf shrub 
species showed a positive correlation between new shoot biomass 
and total ramet biomass (V. myrtillus it 2 = 0.77, P<0.001; V. 
gaultherioides R 2 = 0.7 1, P < 0.00 1 ; E. hermaphroditum R 2 = 0.52, P< 
0.001). 

The number of new shoots per ramet in V. myrtillus was 
significandy greater in warmed plots than in unwarmed plots 
(Fj 9 = 6.466, P= 0.032) and slightly greater under pine than under 
larch (marginally significant, F 1 j S = 3.51, P= 0.077). In addition, 
the mass per individual shoot was larger in warmed plots than in 
unwarmed plots (F ]9 = 9.79, P= 0.012) and slightly larger under 
larch than under pine (marginally significant, Fjjg= 4.14, 
P— 0.060). Therefore, V. myrtillus not only produced more shoots 
in warmed plots, but these individual shoots were longer and 
heavier. Vaccinium gaultherioides produced slighdy more new shoots 
in warmed plots than in unwarmed plots (marginally significant, 
iV >9 =3.74, P= 0.085) and under pine than under larch 
{Fij3= 5.91, P= 0.030), but the mass per individual shoot did 
not differ between warming treatments or plot tree species. 

Growth rings of Vaccinium myrtillus 

We did not find significant differences in ramet age between 
warming treatments (8.4±0.5 years old for unwarmed and 
7.9±0.3 for warmed plots; F y>g = 0.58, P= 0.467; Fig. S2). 
However, ramets were younger under larch (7.4±0.4 years old) 
than under pine (9.0±0.4; F U8 = 10.13, P=0.005). There was a 
marginally significant warming x tree species interaction 
{Fi,is = 3.28, P= 0.087): on average, ramets were older in 
unwarmed plots with pine than in both warmed and unwarmed 
plots with larch. 
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Figure 1. Vegetative and flowering phenology for the three dwarf shrub species studied. Vegetative and flowering phenology for the 
three dwarf shrub species studied. Circles represent the mean day of year (±1 SE) of each phenophase for the four different soil warming (open 
circles, unwarmed; closed circles, warmed) and plot tree species (larch or pine) combinations (n = 10). Crosses show significant differences between 
plot tree species (P<0.05) and marginally significant differences (0.05<P<0.10) between soil warming treatments are given by asterisks in 
parentheses. Colons show treatment interactions. Note that in some cases symbols for different soil warming treatments overlap completely. Dotted 
lines show the date when the warming treatment started in 2012 and arrows on the x-axes show the mean snowmelt date for all plots with the same 
tree species. 
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Repeated measures analyses showed that, averaged across all 
treatment years, V. myrtillus ring width was 4 1 % greater in warmed 
plots than in unwarmed plots (_Fi :9 = 16.45, P= 0.003, Fig. 4). 
Treatment year and the warming x year interaction also had 
significant effects (i^? 97 =4.67, P= 0.001 and F 4:2g7 =3.32, 
P— 0.011, respectively), as well as the pre-treatment ring width 
covariate (^144= 17.26, P<0.00\). Plot tree species and interac- 
tions between tree species and the other fixed effects did not 
significantly influence xylem ring width (P>0.17). Analyses of 
individual years showed that warming had a significant effect on 
ring width in all years (P<0.024) except for the first year of 
treatment (Fig. 4). 

Although xylem latewood formation was not completed when 
we collected the samples in 2012, we found a significant 
correlation between shoot increment length and ring width in 
that year (R =0.157, P<0.00 1), as well as between ring width and 



biomass of new shoots (it =0.088, P<0.001) and biomass of 
leaves (R 2 =0.068, P<0.001). We tested for an age-biomass 
correlation and results corroborated that older ramets were 
heavier (it 2 = 0.24, P<0.001). 

Discussion 

Increased above-ground growth after six years of soil 
warming 

Six years of soil warming led to a sustained growth stimulation 
of V. myrtillus, and the greatest response of new shoot increments 
occurred in the final two years. This increasing response over time 
contrasts with studies where positive responses of V. myrtillus and 
other alpine and Arctic plant species to warming were transient 
and stopped after four or five years [12,13,49]. In addition, our 
findings indicate that other factors that can potentially limit 
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Figure 2. Soil warming effect on dwarf shrub annual shoot 
increment length. Soil warming effect on dwarf shrub annual shoot 
increment length from 2007 until 2012, the entire duration of the soil 
warming experiment. Data through 2009 were presented in Dawes et 
al. (201 la). The warming effect was calculated as the ratio of the mean 
shoot increment length of all warmed plots to the mean of all 
unwarmed plots, pooled across plots containing a larch or pine tree. 
Error bars represent ±1 SE of the ratio. The dashed line shows the 
significant warming effect on V. myrtillus averaged for 2007-2012. The 
asterisk shows significant differences between temperature treatments 
(P<0.05). Pre-warming ratios are shown in the shaded region (2005- 
2006) and the dotted line is drawn through the average of these two 
points, which indicates the mean warmed to unwarmed ratio before 
treatment began. 

doi:1 0.1 371 /journal.pone.01 00577.g002 

growth, such as nutrient or light availability, did not constrain the 
warming response of V. myrtillus in later years of the experiment 
[12]. Earlier studies at the same experimental site found more 
severe freezing damage of V. myrtillus under warming than in 
control plots [36,37], but our results show that these freezing 
effects were not large enough to substantially offset the enhanced 
growth over six years. The other two dwarf shrub species showed 
at least a slight trend of increased shoot increment length with soil 
warming in the last years of the study, suggesting that the response 
of these two species to climate warming may be considerably 
delayed. 

The species-specific warming effects on above-ground biomass 
production in our experiment contrast with a two-year study with 
open-top polythene tents in subarctic Sweden, where all three 
study species in common with our study showed a greater above- 
ground biomass under warming [50] . A five-year study at the same 
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Figure 3. Above-ground biomass partitioning of the three 
dwarf shrub species studied. Above-ground biomass partitioning of 
the study species for each soil warming and plot trees species 
combination (mean values +1 SE, n = 10). Asterisks show significant 
differences between soil warming treatments and crosses show 
significant differences between plot tree species (P<0.05). For 
Empetrum hermaphroditum only, leaves and new shoots are both 
included in "New shoots". The y-axis scale varies across species to 
emphasize differences between treatments. 
doi:1 0.1 371/journal.pone.01 00577.g003 



site in the subarctic by Press et al. [51] similarly showed that 
higher temperatures increased V. uliginosum (comparable with V. 
gaultherioides) biomass. Moreover, Empetrum nigrum shoot production 
and shoot increment length increased under two years [31] and 
five years [32] of warming with open-top chambers in alpine heath 
in Japan. The different methodologies applied (air warming by 
open-top chambers vs. soil warming by heating cables) might 
explain some of the differences between studies. In addition, 
different plant community composition and dynamics might play 
an important role in explaining the different results across study 
sites. 

Vaccinium myrtillus has a wider elevational distribution (colline to 
alpine) than V. gaultherioides and E. hermaphroditum (subalpine to 
alpine) [52]. Its presence in lower elevational areas indicates that 
V. myrtillus might be adapted to higher temperatures than V. 
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Figure 4. Xylem ring widths of Vaccinium myrtillus. Xylem ring 
widths of Vaccinium myrtillus from 2007 to 2011 for warmed and 
unwarmed plots, pooled across plots containing a larch or pine tree 
(mean values ±1 SE, n = 20). Asterisks show significant differences 
between soil warming treatments (P<0.05). Pre-warming values are 
shown in the shaded region (2005-2006). 
doi:10.1371/journal.pone.0100577.g004 



herioides and E. hmnaphroditum [1 1] and thus may be a stronger 
competitor under increased temperatures. In fact, a 22-year 
experiment with open-top greenhouses in the Swedish tundra 
reported no effects on E. hermaphroditum growth under 4°C of air 
warming [13], whereas in the same experiment, Graglia et al. [53] 
found positive effects on this species after six years of 2.5°C 
warming but not 4°C warming. The mean daily air temperature 
increase with open-top chambers is around 1-2°C [54] whereas 
the soil temperature increase by our heating cables is around 4°C, 
which is similar to the air temperature increase reached by 
Campioli et al. [13]. As suggested in their study, E. hermpahroditum 
might suffer from heat stress under high temperatures, responding 
positively only to smaller increases. 

Phenology: effects of growing season temperature, 
snowmelt and plot tree species 

The lack of vegetative phenological responses to soil warming of 
the three dwarf shrub species studied and the strong effect of 
snowmelt date on early-season phenology in 2012 are consistent 
with results from the first three years of the experiment [37]. Some 
plant species have been shown not to respond to spring warming 
experiments [55] because they either do not respond to climate 
warming [55,56], because they are more sensitive to changes in 
late winter temperatures [55] or because they primarily respond to 
other cues such as photoperiod [57]. In our experiment, the lack of 
warming effects on the vegetative phenophases might be partially 
due to the fact that vegetative development, especially for the two 
Vaccinium species, started directly after snowmelt, which coincided 
with the start of the warming treatment. Phenophases occurring 
later in the season (i.e. flowering phenophases in Vaccinium spp.) did 
not depend on the snowmelt date and were slighdy influenced by 
the warming treatment, which is consistent with patterns found for 
tundra dwarf shrub species in subarctic Alaska [58]. 

Our study might have underestimated phenological responses to 
future climate warming because only growing season temperatures 
were altered [55,59], whereas late winter temperatures, which can 
be key determinants of plant phenology [15,55], were not 
experimentally manipulated. Nevertheless, other studies that only 



altered growing season temperatures found that the same species 
as in our study showed phenological advances under air warming 
[18] or that other alpine species showed a lagged response after 
multiple years of warming [10]. The distinct methodologies 
applied in these experiments (e.g. soil vs. air warming), microcli- 
mate conditions, different plant community types and genetic 
variations between populations might explain these differences 
[11,60,61]. 

The tree species present in the plot had a greater effect on V. 
myrtillus phenology than soil warming. The phenological advance- 
ment under larch could be due to a greater light incidence under 
this species at the beginning of the growing season, which could 
act as a phenological cue [62]. 

Increased growth in V. myrtillus decoupled from 
phenology 

The increased growth and biomass production found in V. 
myrtillus with soil warming did not correspond to a phenological 
shift or to the date of snowmelt, which indicates that V. myrtillus 
can produce more biomass under higher growing season 
temperatures even without a springtime extension of the active 
season. The peak-season harvest meant that we could not check if 
V. myrtillus phenology at senescence time was delayed by the 
warming treatment, an effect that has been reported for tundra 
plants in experiments with open-top chambers [63,64]. Delayed 
leaf senescence could have influenced V. myrtillus growth in the late 
season (e.g. secondary growth) or in the following year. However, 
monitoring of leaf senescence in 2008-2009 showed no differences 
between warming treatments in this species (Dawes, unpublished) 
and late-season phenophases of alpine plants are generally 
considered to be more controlled by photoperiod than by 
temperature [65]. 

As advanced phenology could not explain the enhanced growth 
of V. myrtillus under soil warming, potential mechanisms for this 
response include increased rates of photosynthesis [66] and tissue 
formation [65] direcdy caused by higher temperatures. Another 
mechanism could be accelerated decomposition and mineraliza- 
tion of soil organic matter under warmer soil, which can improve 
nutrient availability for plants in alpine and Arctic environments 
where low temperatures tend to constrain these processes [49,66- 
68] . Soil inorganic nitrogen content increased during the first three 
years of our soil warming experiment [1 1], suggesting that indirect 
effects of soil warming via an enhanced nutrient cycling played a 
role in the V. myrtillus growth response. 

Above-ground biomass increase under soil warming 

To our knowledge, this is a pioneer study on assessing the effects 
of experimental soil warming on the age and xylem ring width of 
V. myrtillus. Our results show that rings were wider in warmed plots 
since the second year of treatment, a response that was maintained 
throughout the six experimental years. In addition, the significant 
positive correlation between shoot increment length and early 
xylem growth (xylem latewood formation was not completed) 
indicates that vessel size in the early xylem growing season is 
tightly related to the elongation of new shoots, which receive water 
from newly-formed vessels [69] . 

Moreover, the positive correlation between shoot increment 
length and both early xylem growth and annual shoot biomass 
production, together with the greater biomass in V. myrtillus in 
warmed plots, provide evidence that the previously reported 
increases in the shoot increment length of V. myrtillus with warming 
were not merely a result of a shift in biomass allocation but an 
overall increase in above-ground biomass production. However, 
below-ground biomass was not measured in this study and thus, it 
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is possible that soil warming led to shifts in biomass allocation 
between below-ground and above-ground compartments [70]. 

Although warmer temperatures led to an increased growth of V. 
myrtillus, the age structure of this species was not affected by soil 
warming. However, there was a lower ramet turnover under pine 
(older ramets) than under larch, and ramet main stem biomass was 
also higher under pine. A possible explanation for these differences 
is that pine provides greater protection against freezing conditions 
at the beginning of the season before needles are developed on 
deciduous larch. This effect would be similar to the facilitation 
exerted by shrubs on young trees [71,72], leading to lower 
mortality rates. Moreover, Dawes et al. [1 1] found less canopy 
shading under pine than under larch, thus the lower main stem 
biomass under larch may be a consequence of lower light 
availability. 

Concluding remarks 

In summary, we found increased growth of V. myrtillus under soil 
warming, a response that was sustained, and even became stronger 
in the case of shoot increment length, over six years of warming. 
The application of dendrochronological techniques showed that 
this increased growth reflected an overall increase in above-ground 
biomass production. Moreover, the lack of an advanced phenology 
of V. myrtillus under soil warming indicated that an extended 
growing period was not necessary for the observed growth 
response. Our results suggest that V. myrtillus will experience a 
larger and more rapid growth benefit from a warming climate 
than V. gaultherioides or E. hermaphroditum and could therefore 
become increasingly dominant in high-elevation treeline environ- 
ments. 
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